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Context. Current atmospheric and evolutionary models for low-mass stars and brown dwarfs rely on approximate assumptions on the 
physics of the stellar structure and the atmospheric radiative transfer. This leads to biased theoretical predictions on the photospheric 
Spectral Energy Distributions of these system, especially when applied to low surface gravity objects such as Pre-Main Sequence 
(PMS) stars, and affects the derivation of stellar parameters from photometric data. 
Aims. Our main goal is to correct the biases present in the theoretical predictions for the near-IR photometry of low-mass PMS 
stars. Using empirical intrinsic IR colors, we assess the accuracy of current synthetic spectral libraries and evolutionary models. We 
investigate how the uncertainty in the intrinsic colors associated with different PMS models affect the derivation of the Initial Mass 
Function of young clusters from near-IR photometry. 

Methods. We consider a sample of ~300 PMS stars in the Orion Nebula Cluster (age^l Myr) with well measured luminosities, 
temperatures and photospheric JHKs photometry. This sample is used as a benchmark for testing both atmospheric and evolutionary 
theoretical models. 

Results. By analyzing the photospheric colors of our sample of young stars, we find that the synthetic JHKs photometry provided 
by theoretical spectral templates for late spectral types (>K6) are accurate at the level of ~0.2 mag, while colors are accurate at 
<0. 1 mag. We tabulate the intrinsic photospheric colors, appropriate for the Orion Nebula Cluster, in the range K6-M8.5. They can be 
OO ' conveniently used as templates for the intrinsic colors of other young (age<5 Myr) stellar clusters. 

Conclusions. The theoretically-predicted JHKs magnitudes of young late type stars do not accurately reproduce the intrinsic ones of 
the Orion Nebula Cluster members. An empirical correction of the atmospheric templates can fix the discrepancies between expected 
and observed colors. Still, other biases in the evolutionary models prevent a more robust comparison between observations and 
theoretical absolute magnitudes. In particular, PMS evolutionary models seem to consistently underestimate the intrinsic near-infrared 
flux at the very late spectral types, and this may introduce spurious features in the low-mass end of the photometrically-determined 
Initial Mass Function of young clusters. 
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ABSTRACT 



1. Introduction to t he substellar regime with the COND and DUSTY mod- 
els dChabrier et all l2000h : lYi et all d2003l) : and iTognelli et ail 

In the last two decades several generations of evolution- (l201ll) . However, considerable differences still exist as the var- 

ary models of pre-main-sequence (PMS) objects have been ious sets of evolutionary tracks show systematic differences 

published with continuous improvements on the treatment in the predicted masses and ages of stars in the H R diagram 

of the hydrostatic structure, radiation/heat transfer and fher- dHillenbrand & Whitdl2004t iHillenbrand et al.ll2008l) . 
modynamic equilibrium. There are at least 7 different fam- The comparison of theoretical predictions with observations 

ilies of published PMS evolutionary calculations that have on the HR diagram requires converting observed quantities, e.g. 

been widely circulated in machine-rea dable formats and tha t spectral type and some photometry in at least two passbands, 

span a suitable range of stell ar masses: ISwenson et alJ (Q994); into physical param eters like effective temperature and absolute 

iD'Antona & Mazzitellil (fT997l) w it h 1998 elect r onic-o nly update luminosity (see e.g. lHillenbrandl[l997l) . This conversion is nor- 

(DM98): iPalla & Stahlerf (fl999h: ISiess et al.1 d2000h (SDFOO); mally done using an empirical, or semi-empirical, calibration. 

lBaraffeetal.l (119981) (BCAH98) with the subsequent extensions It is also possible to take the complementary approach of de- 
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riving observational data from theoretical predictions. For ex- 
ample, plotting the photometry of a young stellar cluster in a 
color-magnitude (CMD) or two-color diagram (2CD) may show 
structures (like e.g. the cluster isochrone) that can be compared 
with the results of evolutionary and synthetic photometry calcu- 
lations. Varying the model parameters one can directly visual- 
ize the effects of metallicity, effective tempe rature scale, surface 
gravity, reddening law, accretion (see e.g. iDa Rio et al.1 BOlOi 
D10 hereafter). 

No approach is without drawbacks. On the one hand, the em- 
pirical corrections needed to derive the physical parameters of 
PMS stars, like e.g. colors and bolometric corrections, may have 
been derived on samples of stars which may not be fully rep- 
resentative of the sample under examination. There is a stand- 
ing tradition, for example, of using for P MS stars the intrin- 
sic co lors of Main Sequence stars (e.g. iKenvon & Hart mann 
1 1995b . The long-known difference between the spectral type vs. 
temperature relations for dwarfs and giants represents another 
source of uncertainty for PMS stars, that many authors overcome 
(or mitigate) th is uncertainty using the "intermediate" scale of 
Luhman ( 1999), which has been found ad equate for PMS stars 
in many studies (e.g., IDa Rio et al.ll2Q10l) . On the other hand, 
the systematic uncertainties of evolutionary models may com- 
bine with those of the atmospheric templates and provide er- 
roneous results. The selection of adequate model atmospheres 
for PMS stars is especially critical for late type stars (M-type 
and later), which rep resent the peak of the Initial Mass Function 
(Bast ian et al.l 1201 CH) and whose spectral energy distributions 
(SEDs) are dominated by broad molecular absorption features. 
Cool atmospheres (T^ <3000 K) host a variety of molecules and 
provide an environment for dust condensation, whose role in the 
heat and radiation transfer through the photosphere cannot be ne- 
glected (lAllardet al.ll200ll 12010). Moreover, for sub-stellar ob- 
jects, convection may involve the photosphere, which therefore 
can no longer be treated as a system in radiative equilibrium. It 
is therefore crucial, particularly for cool atmospheres, to include 
in the synthesis codes the largest number of molecular lines in 
order to accurately reproduce the radiation transfer through the 
atmosphere. 

These difficulties are aggravated by observational limita- 
tions. It is hard to obtain high quality data for a statistically 
significant sample of PMS stars of comparable age and dis- 
tance. Even in the solar vicinity, rich and young stellar clus- 
ters are typically affected by large and inhomogeneous redden- 
ing, being still ensh rouded in their parental molecular cloud 
(lLada&Ladall200l . while the stars often show accretion ex- 
cess. Older systems tend to be spatially spread and affected by 
membership uncertainties. In this respect, due to its relat ively 
low foreground extinctio n (Ay < 3,IScandariato et al]|201 lb and 
vicinity (d = 414 pc, iMenten et all 120071) . the Orion Nebula 
Cluster (ONC) provides a unique opportunity to analyze the in- 
trinsic colors of its members and assess the accuracy of PMS 
models. 

In the optical wavelength range, a recent attempt t o calibrate 
empir ically the colors of PMS stars is presented in lDa Rio et al.1 
(2009). In that work, based on optical BVI photometry of the 
ONC, the authors show how present families of synthetic spec- 
tra fail in matching the observed colors, and present a correction 
based on their data. In D10, they further refine this calibration, 
limited to the / band and 2 medium bands at A ~7700 A, by ex- 
plicitly calibrating colors as a function of T e s, and by decreasing 
the lowest T e ff limit down to ~2800 K. 

The goal of this paper is to use the extensive set of spectral 
types and photometric data available for the ONC to test the the- 



oretical models and to empirically determine the intrinsic (pho- 
tospheric) JHKs magnitudes and colors of the cluster members 
as functions of T e g . The intrinsic NIR colors of PMS stars in 
Orion we derive will be also appropriate for the (possibly ideal) 
cluster isochrone, and applicable to young (age<5 Myr) systems 
in general. In Sect. [2] we present the selection of our sample of 
stars, based on the latest spectral characterization and our re- 
cently published NIR photometry. In Sect. [3] we refine our list 
to select the subset of stars most suitable for our purposes. By 
means of this samp le, in Sect . H we test the most recent atmo- 
spheric model of Allar d et alj ( 20101) . and in Sect. [5] we derive 
the average colors of the cluster. Finally, in Sect. [6] we discuss 
our result and we compare them to the current theoretical pre- 
dictions. 



2. The data set 

W e base our anal y sis o n the ONC NIR photometric catalog 
of iRobberto et alJ (120101) . This catalog contains JHKs pho- 
tometry in the 2MASS system for ~6500 point-like sources, 
spread over an area of ~30'x40' roughly centered on fl'Ori-C 
(RA=05''35 ffl 16.46 s , DEC=-05°23'23.2"), down to JHK S ~18. 
W e cross-match the N IR catalog with the list of BVI photometry 
of|D a Rio et al] d2009l) . These NIR and optical catalogs are prod- 
ucts of the HST Treasury Program on the Orion Nebula Cluster 
(HST GO-10246). They originate from simultaneous observa- 
tions of the cluster, and therefore are largely immune from the 
uncertainties introduced by source variability when multi-band 
photometry is collected at different epochs. 

IDa Rio et al.l (120121 D12 hereafter) provide photospheric lu- 
minosity, interstellar extinction and accretion luminosity for 
~ 1,200 stars, whose spectral types are known e ither from pre- 
vious spectroscopic surveys (iHillenbrandll 19971) or from obser- 
vations of the narrow -band photometric index of TiO, either new 
or from D10. 

In order to extend our analysis down to later M-subtypes 
and probing the Brown Dwarfs (BDs) spectral ra nge, we com- 
pleme nt our list with the spectral classification of iRiddick et al.l 
(2007, R07 hereafter) for 45 M-type stars not in our catalog. 
The authors provide extinctions, luminosities and spectral types, 
together with JHKs photometry from their previous NIR photo- 
metric survey dLucas et alj|2005l) . 

We also add 51 m ore stars from the spectroscopic study of 
ISlesnick et al.l d2004l S04 hereafter), which also provides extinc- 
tion, luminosity and spectral type of the s tars. H Ks photome- 
try is taken from Hillenb rand & Carpenter! (|2000), whereas no J 
photometry is available. 

In the last two cases, individual extinctions are derived spec- 
troscopically, while the photometric system of the observed NIR 
photometry is consistent to the 2MASS system at a level of 
<0.1 mag. This uncertainty is lower than the scatter introduced 
by stellar variability , which typically pro duces fluctuations of the 
order of >0.2 mag dCarpenter et al .11200 ll) . 

We correct the observed NIR p hotometry for interst ellar ex- 
tinction using the extinction law of lCardelli et al.l dl989l) and the 
values of Ay for individual stars derived spectro-photometrically 
by D12, R07 and S04. Spectr al types are converted into effec- 
tive temperature r e ff using the iLuhman et aQ (120031) conversion 
scale, commonly adopted for PMS stars (Table Our list thus 
contains 1321 stars with spectral type, luminosity L, extinction 
Ay and extinction-corrected 2MASS-calibrated NIR photome- 
try; 92% of the sample comes from D10, the remaining part 
coming from either R07 or S04. 
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3. Extraction of the reference stars 

In order to identify a "master" sample of PMS stars that may 
more adequately represent the intrinsic colors of the ONC pho- 
tospheres, we need to clean our sample from field stars, outliers 
and members with NIR fluxes contaminated by close-in com- 
panions, circumstellar emission or strong activity. 

3. 1 . Age spread 

First, we reject outliers in the HR diagram. Using the list of 
r e ff and L compiled in the previous section, we derive ages and 
masses through a comparison of their positions in the HR dia- 
gram with the evolutionary models of DM98 (Fig. [TJ. The ad- 
vantage of the DM98 models with respect to the others is that 
they nicely cover the mass and age ranges spanned by the ONC, 
while, e.g., the BCAH98 models are available for age>l Myr, 
consistently older than the youngest stars in the HR diagram 
shown in Fig.Q] 

The distribution of derived logarithmic ages, plotted in 
Fig. [2] indicates that according to this model the characteristic 
age of the ONC members is ~l-2 Myr. Aiming at a large sta- 
tistical sample, representative of the whole cluster, we conserva- 
tively retain all stars with age within the full width at one tenth 
of the maximum (4.6 < Log(age) < 7.2), rejecting all those that 
appear too young or too old with respect to the average popula- 
tion. 

This selection thus rejects stars too bright compared to the 
average cluster, i.e. stars in the foreground, intrinsically too 
young or possibly associated to circumstellar excess or close 
companions. 

This selection also discards the stars too faint compared to 
the average cluster, li kely belonging to o lder stellar associations 
in the Orion complex (Brown et al .il994h or to the galactic back- 
ground population dDa Rio et al.ll2012l) . 

Using other models would have provided different values for 
individual ages, mean cluster age and lower/upper age limits, but 
our criteria would have basically selected the same sample. 

3.2. Stellar multiplicity 

Tobin et al] d2009l) find that the frequ ency of spect roscopic bi- 
naries in the ONC is -10%, while iKohler et al.1 d2006l) and 
Reipurth et al. (2007) indicate that the frequency of visual bina- 
ries in the cluster is < 10% at separations >60 AU. In particular, 
Kohl er et al.l ((2006) find that the frequency decreases down to 
~5% across the whole cluster for low-mass systems. 

To clean our sample from unresolved multiple systems, we 
checked the HST/ACS image s taken for the Orion Treasury 
Program (Rob berto et al.ll2012l) . This check allows us to discard 
40 stars with projected companions at distances >50 mas (the 
angular resolution of the images), unresolved in the ISPI images. 

For angular distances <50 mas, corresponding to ~20 AU 
at the distance of the ONC, any multiple systems would remain 
unresolved. We thus expect that the fraction of unresolved mul- 
tiple low-mass sys tems in our sample is <10%, as indicated by 
Tobi n et alJ 42009), and does not influence significantly our re- 
sults. 

3.3. Extinction 

The extinction estimates in our master sample are derived ei- 
ther photometrically (D10) or spectroscopically (S04, R07). In 
both cases, the observations are compared to either theoretical 



i i i i i 




7.5: 

' ' ' ' ' - 



3.9 3.8 3.7 3.6 3.5 3.4 3.3 
Log(T eff ) (K) 

Fig. 1. HR diagram of our list of stars superimposed to the evo- 
lutionary model computed by DM98. Solid lines show a few 
isochrones at ages included in our selected age range (the log- 
arithmic age in year is indicated), while dashed lines show the 
tracks at fixed mass. Dots, crosses and squares represent the stars 
from D10, R07 and S04 respectively. The master sample ana- 
lyzed in Sect.|5]is shown in red. 

1 .o ; : 




Log(age) (yr) 

Fig. 2. Age distribution of the stars based on the DM98 evolu- 
tionary models. This plot indicates that the ONC is ~l-2 Myr 
old. The two dotted lines outline the age range spanned by the 
analyzed sample of stars. 



or empirical templates. This approach is somehow inaccurate as 
the effects of surface gravity becomes relevant at spectral types 
later then ~M3. In order to avoid these uncertainties, we re- 
tain the ONC members poorly embedded in the nebula, selecting 
the stars in our master sample characterized by low extinctions 
(A v <0.5). 

3.4. Circumstellar Activity 

As shown bv lCarpenter et alJ (l2002l) . accreting stars are charac- 
terized by variability up to 0.2 mags in the NIR. Thus, to avoid 
the uncertainties introduced by accretion, we reject those stars 
with significant circumstellar activity. 
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Fig. 3. Distribution of the accretion luminosities / in our selected 
sample. Values / <-6 must be regarded as null detection of on- 
going accretion, and are kept in the analyzed sample. On the 
contrary, stars with / >-6 show accretion evidences and are dis- 
carded, in order to avoid any contamination to the extinction- 
corrected photometry. 



For the subsamp le of stars with BVI photometry taken from 
iDa Rio et alj d2009t) . D10 provide indications on ongoing accre- 
tion through the parameter I defined as: 



l = Log\ 



where L tot is the total luminosity of the system and L accr is the 
accretion luminosity. We discard stars with I >-6 (Fig. [3]), in 
order to isolate only those showing no evidence of accretion. 

We also exclude the stars in the R07 and S04 sublists with 
spectroscopic evidence of disk activity, mainly Call emission 
lines which are commonly attributed to the presence of disk 
winds. 

3.5. NIR excess from the inner disk 

Circumstellar disks aroun d young PMS star s are known to ra- 
diate in the NIR d oma in dMever et alJll997t ICieza et alJ I2005L 
iFischer et al.ll2.Ql il andlHillenbra nd et al.l (I1998) for the specific 
case of the ONC), contaminating the NIR continuum emission 
of the central stars. To reject sources with photometry contam- 
inated by strong disk excess, we plot in Fig. [4] the NIR 2CD of 
our 395 stars, comparing their extinction-corrected colors to the 
1 Myr isochrone of DM98. 

Whereas the majority of stars appear scattered around the 
DM98 isochrone, the redder sources (top-right corner of the plot) 
appear systematically shifted in the region corresp onding to the 
Classi cal T-Tauri Stars (CTTS) locus proposed by iMever et aD 
(Il997l) . These dereddened colors cannot be accounted for by ei- 
ther inaccuracies in the theoretical 2CD or by measurement er- 
rors, and must therefore be attributed to circumstellar disk emis- 
sion. Taking the theoretical isochrone as a proxy for the photo- 
spheric colors of the selected stars, we discard the stars with 
extinction-corrected photometry deviating from the isochrone 
for more than twice their photometric uncertainties, and we keep 
the other stars, whose infrared photometry is consistent with the 
theoretical models within the uncertainties. 



Fig. 4. 2CD of our sample of stars, compared to the 1 Myr old 
model provided by DM98 and computed in the 2MASS system 
(pink solid line). The observed colors do not strictly follow the 
theoretical model, indeed the y occupy a r e gion e longated along 
the CTTS locus provided by Mey er et aD ([1997) (green dashed 
line). Cir cles and crosses repr esent stars with spectral types pro- 
vided by IDa Rio et al.l (120121) and R07 respectively. Gray error 
bars represent stars with observed colors not compatible with 
the model at the 2<x level: they are rejected in our analysis. 



To check the robustness of this selection criterion, we cross- 



match our full NI R catalog dRobberto et all 12 010) with the 
YSOVAR catalog dMorales-Calderon et alJ 1201 lh . containing 
Spitzer/IRAC photometry of the ONC stars. Based on the ob- 
served mid-infrared colors, this catalog also classifies the stars 
in Class II stars (i.e. PMS stars with circumstellar disks) and 
Class III stars (i.e. PMS stars with debris disks, weakly radiating 
in the NIR). 

We analyze the cross-matched catalog by means of the NIR 
2CD, as shown above, and we find that >60% of Class II stars 
are consistently displaced in the diagram along the CTTS locus, 
while >70% of Class III stars have NIR colors compatible with 
the theoretical isochrone. This indicates that the majority of the 
selected stars are Class III stars, while the Class II stars in the 
retained sample likely have negligible NIR excess, compared to 
the photometric uncertainties. 

This is the last step of our selection process. We end up 
with a master sample made of 150 bona-fide cluster members 
with 2500 K<r e ff<5000 K, poorly extincted and without sub- 
stantial evidence of non-photospheric activity, either accretion 
or circumstellar disk emission. 



4. Accuracy of the lAllard et al.ll|2010|) atmospheric 
model 

In this section, we take advantage of the list of stars assembled in 
Sect. [2] to investigate t he accuracy of the s ynthetic atmospheric 
templates provided by lAllard et aH (l2010h . which nicely cover 
the parameters space spanned by our sample of stars. 

The assignment of spectral types of late-type young stars is 
usually performed in the optical-red domain, as this wavelength 
range provides the best-studied sets of spectral lines for clas- 
sification and is least affected by veiling from accretion and/or 
circumstellar emission. With spectral types at hand, it is then 
possible to obtain a prediction of the photospheric fluxes in a 
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certain passband by convolving the corresponding filter profile 
with an appropriate template spectrum, once 7^, the logarith- 
mic surface gravity log g and the metallicity are known. 

For each star in our master sample we already have the ef- 
fective temperature at hand, and we assume t he standard so- 
lar me tallic ity for the ONC, as indicated by, e.g.. lD'Orazi et al.l 
(2009) and Bi azzo et al.l (1201 ll) . In order to compute synthetic 
photometry, we thus need an estimate of log g. Having the lu- 
minosity L and Teg, we first compute the stellar radius R* us- 
ing the Stefan-Boltzmann's law L = 4jtRIct s B T^ S . We also de- 
rive stellar masses M* placing the stars in the HR diagram and 
comparing their position with the theoretical model of DM98 
(Fig-ID - We refine the mass estima tes taking into account the re- 
sults of iHillenbrand & White! d2004l) . who analyzed the discrep- 
ancies between track-predicted masses and dynamical masses 
for a sample of PMS stars in binary systems. They found that 
the DM98 models tend to underestimate masses by ~20%. We 
then increase by 30% our track-predicted mass estimates, and 
we use the resulting values to compute the surface gravities g 
following the relation g - GM*jR\. 

For each star, we then compute synthet ic 2MASS N I R pho - 
tometry using the atmospheric models of All ard et all (|2010), 
with the appropriate values for the effective temperature and 
sur face gravity, and t he distance d=414+7 pqj as derived 
by iMenten et alJ (|2007). The synthetic magnitudes and colors 
^synthetic are then compared to the corresponding extinction- 
corrected values X ex ,j„ c ,j „- correcte d, computed taking into account 
the Ay estimates from optical spectro-photometry (see Sect. [2]). 
The differences 



AX(T eff ) = X, 



extinct ion -cor reeled 



synthetic 



(T eS ) 



(1) 



represent our proposed corrections to the synthetic NIR colors. 
These are values to be added to the synthetic colors to derive the 
correct intrinsic colors of the stars. 

The distribution of stars in terms of T e ff allows to perform 
a smooth non-parametric fit of the AX from r e ff~2500 K up to 
T e ff ~4200 K, while the number of stars at higher temperatures 
is too low to allow for a robust fit (see Fig. [TJ. The fitted AH, 
A(J-H) and A(H-Ks ) for the entire master sample are shown in 
Fig.|5]and reported in TableQ] as a function of the effective tem- 
perature. They indicate that the H magnitude is underestimated 
by the model (i.e. the model systematically provides brighter H 
magnitudes) by <0.2 magnitudes, with an uncertainty <0.1 mag- 
nitude. The same is true for the colors, with a smaller offset of 
the order of <0. 1 magnitudes, in the red and blue directions for 
the J-H and H-K$ colors respectively. 

One possible explanation for these discrepancies is that the 
atmospheric model provides H fluxes larger than the observed 
ones, while providing fairly good predictions on the JK$ fluxes. 
From an observational point of view, this hypothesis is consis- 
tent with Fig. |6l where we show that the atmospheric model 
generally provides consistent predictions on the J-K$ colors. In 
order to address this issue, in Fig. [7] we compare the synthetic 
M7 spectrum of Allard et al. (2010) with the spectra of two M7 
stars observed bvlMuench et alJ d2007l) in IC 348 (age~2 Myr, 
iLuhman et al.ll2003l) . This plot confirms that, despite the im- 
provements in the latest release, the model still does not fully 
account for the continuum opacity of water vapor, which signifi- 
cantly affects the NIR SED, and in particular the H-band contin- 
uum, of cool stars. 
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1 The uncertainty on the distance of the ONC negligibly propagates 
onto the uncertainties on the observed photometry. 



Fig. 5. AH, A(J - H) and A(H - K s ) corrections to the syn- 
thetic colors (from top to bottom respectively). Circles, crosses 
and squares repres ent stars with spectral types provided by 
iDa Rio et all (120121) . R07 and S04 respectively, while the error 
bars take into account the uncertainties in both the extinction- 
corrected photometry and the synthetic photometry. Red solid 
lines show our smooth outlier-resistant fits, together with the cor- 
responding 95% confidence bands (dashed red lines). 



Another possibility is that the atmospheric models do not ac- 
count for the collision induced absorption (CIA) from molecular 
hydrog en H2. Theoretical models computed by iBorvsow et al.l 
dl997l) show that CIA is relevant for the infrared spectrum of 
cool subgiants: radiative energy is absorbed in the infrared and 
re-emitted at visible wavel engths, leading to the typical trian- 
gular shape of the H-band dAllers et alJl2007b iKirkpatrick et all 
2008, see also Fig.Qand Fig.|9]in this paper). 



4.1. Systematic uncertainties 

The corrections we have derived are based on the reddening law 
used to correct the observed photometry for extinction and on the 
accuracy of the mass estimates. The results may be influenced by 
the assumed ratio of total to selective extinction Ry (Sect. l4~T7Tb 
and by the accuracy of the mass (and surface gravity) estimates 
derived from the DM98 theoretical model (Sect. 14. 1.21 ). 



4.1 .1 . Effects of the assumed R v 

In our derivation of the intrinsic NIR colors of stars in the ONC, 
we have assumed that the total to selective extinction is Ry = 
Ay/E(B - V)=3.1, which is the typical Galactic value. However, 
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Fig. 7. Synthetic spectra of a M7 PMS star (solid line) as p rovided bv lAllard et alj (l2010l) . together with the observed spectra of two 
young M7 stars (in dots and dashes, iMuench et al.l l2007). The blue, green and red lines stand for the transmission curves of the 
2MASS JHKs filters respectively. We find good agreement between the model and the observations in the JK$ bands, while in the 
fi band the model overestimates the radiated flux (or underestimates the opacity). 



Table 1. Corrections to synthetic colors for LMSs and BDs in 
the ONC. 
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Notes. (<!) Temperature scale of Luh man et al] (12003). 
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Fig. 6. A(J — Kg) corrections to the synthetic colors. Symbols 
and colors are the same as in Fig. [5] 



it has been proposed that the reddening law toward the OB stars 
of the Orion association can b e better parametrized by Ry ^ 5.5, 
typical of larger dust grains ( Johnson! 1967t Costero & Peimbert 
19701; iBaldwin et al.l fl99ll lOsterbrock et all 1 1992b iGreve et a" 
1994HBlagrave et alj|2007l) . 

We thus investigate how our results vary with increasing Ry 
considering the subsample of stars of D10 (corresponding to 



~92% of our collected catalog down to spectral type M7), for 
which the authors provide the stellar parameters derived assum- 
ing Ry=5.5. We do not include the list of stars of R07 and S04 
in this analysis as the authors do not derive the stellar parame- 
ters assuming Ry =5.5. By consequence, we can derive the AX 
corrections only down to the M7 spectral type. 

The new set is analyzed following the same strategy de- 
scribed in Sect. [3] Using the new set of the stellar parameters L, 
Ay and I, we compile a catalog of 1 17 stars and derive new cor- 
rections to the NIR synthetic colors (TableQ]). In Fig. [8] we com- 
pare these new corrections to the ones obtained with Ry=3.l, 
finding that the two sets of corrections are consistent within the 
errors. 

Moreover, as shown by IScandariato et alj d201 lh . the 
foreground Orion Nebula generally provides magnitudes 
Ay <2 mag. Thus, the cut at low extinctions discussed in 
Sect. l3.3l excludes all those stars deeply embedded in the nebula. 
We also remark that the same Ry holds for higher extinctions 
within the Orion Nebula as in dicated by e.g. Indeb etouw et alj 
(120051) and lDa Rio et alj d201dh . 

For these reasons, and considering the weak dependence of 
IR extinction on different values of Ry, from now on we will 
adopt for simplicity the case Ry-3.1. 

4.1.2. Surface gravity estimate 

Our analysis of the synthetic photometry has been based on a 
particular family of evolutionary models, DM98, used to de- 
rive the stellar mass and therefore the surface gravity. This as- 
sumption may affect especially the NIR continuu m of M-type 
stars, shaped by the water vapor opacity profile dAllard et al] 
2010), which is sensitive to surface gravity. In order to reduce 
the uncertainties in log g related to the assumption of a partic- 
ular model, in Sect. [4] we increas ed the DM98 track-predicte d 
masses by ~30%, as indicated by Hillenbrand & White! (|2004). 
Since mass and surface gravity are linearly correlated, this cor- 
rection increases \ogg by ~0.16 dex with respect to the value 
derived using the DM98 model masses at face value. 

To further investigate the effects of surface gravity in our 
analysis, in Fig.|9]we compare the synthetic near-IR spectra of an 
Ml and M7 star both 1 Myr old according to DM98. The dashed 
and solid lines show the spectra before and after the mass correc- 
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Fig. 8. Comparison of the AH, A(7 - H) and A(H - K s ) correc- 
tions (from top to bottom respectively) derived assuming Ry-3.1 
(blue lines) ans Ry-5.5 (red lines). The dashed lines represent 
the 95% confidence bands. 



tion, respectively. For the M7 star the correction has no impact. 
On the other hand, the HKs continuum of the Ml star decreases 
by <0.01 dex, while no significant changes are present in the / 
band. In terms of synthetic photometry, the increase of surface 
gravity leads to fainter HKs magnitude and bluer J-H colors by 
<0.03 magnitudes, which is less than the uncertainties on the at- 
mospheric model (see, e.g., Fig. O. Thus, fine tuning the \ogg 
values has a negligible effect on the synthetic photometry and, 
by consequence, the differences between predicted and observed 
photospheric colors are not dominated by the uncertainties in the 
evolutionary models, as long as the track-predicted masses devi- 
ate by <30%. 

If using different recipes for deriving the stellar mass, and 
therefore surface gravity, has a negligible effect on the synthetic 
spectra of PMS stars, the discrepancies we found between model 
predictions and extinction-corrected photometry (Table [1} must 
be largely attribute d to the re s idual uncertainties of the atmo- 
spheric models of lAllard et alj d2010l) in the NIR. 



5. The empirical NIR isochrone of the ONC 

By taking advantage of the sample of ONC stars compiled 
in Sect. [2] we can empirically derive a representative NIR 
isochrone of the cluster, assuming Ry=3.1. In Fig.[l0]we show 
the extinction-corrected H, J-H and H-Ks colors of the sam- 
ple of stars analyzed in Sect. [3] as functions of the effective 



Fig. 10. Intrinsic H, J-H and H-K s colors of stars in the ONC 
as function of T e ff (from top to bottom respectively). Circles, 
cr osses and squares rep resent stars with spectral types provided 
bv lDaRio etalj d2012l) . R07 and S04 respectively, while the er- 
ror bars take into account the uncertainties in both the extinction- 
corrected photometry and the synthetic photometry. Red solid 
lines show our fits, together with the corresponding 95% confi- 
dence bands (dashed red lines). 



Table 2. Empiric intrinsic colors for LMSs and BDs in the ONC. 



SpT 


7 a 
Jeff 
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J-H 


H-K s 


M8.5 


2555 


15.21 


0.52 
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M8 


2710 
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0.53 


0.45 


M7 
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M5 
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M4 
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M3 
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0.24 


M2 
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0.20 


M0 


3850 


10.96 
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0.17 


K6 
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10.56 
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Notes. <D) Temperature scale of Luh man et alJ (T2003). 



temperature. Following the same strategy adopted in Sect. [4] 
we smoothly fit the intrinsic colors from r e jf~2500 K up to 
J e ff~4200K (3.40<Log(r e jf)<3.62) using a local polynomial re- 
gression smoother. The derived intrinsic NIR colors of the ONC 
are given in Table [2] 
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Fig. 9. Synthetic spectra for a Ml and a M7 star (top and bottom spectrum respectively) as provided by lAllard et alj d2Q10t) . For 
each spectral type, the dashed line shows the spectrum computed with logg provided by the 1 Myr old model of DM98 (3.63 and 
3.33 respectively), while the solid line shows the same spectra with logg artificially increased by 0.16 dex (for the M7 star the two 
spectra overlap). The blue, green and red lines stand for the transmission curves of the 2MASS JHK$ filters respectively. We find 
that the M7 spectrum weakly depends on log g. On the contrary, for the Ml star the continuum spectrum in the HK$ bands decreases 
with increasing log g. No changes are seen in the J band. 



The magnitudes and colors presented in Table [2] represent 
our best estimate of the locus in the NIR CMDs along which the 
ONC stars should be located if one removes the effects of extinc- 
tion and NIR excess, and there is no intrinsic color dispersion 
among stars. We underline that this is an average, ideal locus. 
The ONC population is spread across the HR diagram (see e.g. 
Fig- [!]) by about ~1.5 dex in luminosity, partially due to effects 
like stellar variability and unresolved companions. Moreover, re- 
cent studies indicate that the l uminosity spread may be ascribed 
to an intrinsic age spread (Hillenbrand 2009; Reggiani et al. 
1201 lb . The observed luminosity spread, ~3.5 magnitude, is re- 
flected in the top panel of Fig. [10] by a spread along the ver- 
tical axis larger than the error bars. Our empirical H magni- 
tudes, therefore, may only be appropriate for an average ONC 
isochrone. 

On the other hand, we have seen that the differences in lumi- 
nosity, while leading to differences in surface gravity, have small 
influence on the photospheric colors (Sect. l4.lT2l ). For this reason 
the average NIR colors in Fig. [10] are better constrained than the 
absolute magnitudes. One can therefore assume that our empiri- 
cal, average colors properly describe those of the ONC sources, 
and of their corresponding isochronal sequences. 

In Fig.QT]we compare our empirically-derived colors for the 
ONC stars with the correspo nding observed average colors of 
giants and dwarfs provided bv lBessell & Bred (Qj>88). This plot 
shows that the J-H and J-K$ colors generally get bluer with in- 
creasing gravity, while the H-K$ color is less affected. This is 
consistent with Fig. [9] where we show that the H and K$ fluxes 
tend to increase with log g, while the / flux is not sensitive to 
surface gravity. Moreover, we remark that the ONC stars, be- 
ing PMS stars, have intermediate gravities between giants and 
dwarfs, and this is reflected by the fact that they have NIR colors 
intermediate between the giants' and dwarfs' colors. 

5.1. Test of methodology 

As a sanity check, we now test our methodology against a sam- 
ple of main sequence stars in the range 5,000 K<r e ff<6,000 K 
(corresponding to the K1-G0 spectral type in the temperature 
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(green squares) colors of the ONC (solid line ), giants (dotted 
lines) and dwarfs (dashed lines), as provided bv lBessell & Brettl 
(1988). 
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Fig. 12. Observed H, J-H and H-K$ colors of main sequence 
stars in the in the solar neighborhood as function of T e ff (from 
top to bottom respectively). Red solid lines show our fits, to- 
gether with the corresponding 95% confidence bands (dashed 
red lines), while the cyan lines represent the synthetic colors 
computed as described in the text. 



scale of Schmidt-Kaler (Il982ln . where theoretical models are 
known to produce robust predictions on p hotospheric colors. W e 
select 77 stars in the PASTEL catalog (ISoubiran et alj 120101) . 
with solar metallicities, surface gravities and temperatures de- 
rived spectroscopically, and JHK$ magnitudes extracted from 
the 2MASS catalog. These stars are located within 30 pc from 
the Sun and, by consequence, the observed NIR photometry is 
negligibly affected by interstellar extinction. 

In Fig. [12] we plot our best fits of the photospheric colors 
of the selected stars as functions of T e ff, obtained as described 
in Sect. [4] We also compare our empirical fits with the syn- 
th etic photometry comput ed combining the evolutionary model 
oflTTA ntona & Mazziteflj (Il998l) with the atmospheric model of 
Allar d et all Q2010). These models have been extensively tested 
against main sequence stars, and have proved to consistently re- 
produce the observed NIR magnitudes and colors for the an- 
alyzed spectral range within the photometric uncertainties. A 
posteriori, we fix the inaccuracies of the atmospheric model de- 
scribed in Sect. [4] adding 0.05 magnitudes to the synthetic H 
magnitudes, thus extrapolating Table Q] towards earlier spectral 
types. 

Figure [12] shows that our algorithm recovers the expected 
magnitudes and colors within the 95% confidence bands. This 
test thus confirms the consistency of our analysis, providing ad- 
ditional confidence to the empirical colors listed in Table [2] 



6. Comparison with theoretical models 

6.1. Magnitudes and colors vs. T eS 

As anticipated in Sect. Q] current evolutionary models provide 
different predictions for the position of the PMS isochrones and 
evolutionary tracks in t he HR diagram, especia lly at sub-solar 
masses and young ages (Iffillenbrand et al.l l2008). This is a con- 
sequence of the difficulties in adequately modeling the convec- 
tive equilibrium in the interior of fast rotating sub-solar young 
stars. 

Alongside with the differences between the adopted input pa- 
rameters and theoretical treatment, theoretical models also differ 
in setting the t — instant of the PMS evolution, i.e. there is 
no absolute calibration of the "clock" counting the age of young 
stars. For this reason, theoretical models may also mismatch. 

The inconsistencies between the various sets of models gen- 
erally persist when they are converted into magnitudes and col- 
ors, for a direct comparison with the data in the CMDs (see 
Fig. fT4T> . As discussed in Sect. I4.1.2I the synthetic magnitudes 
depend weakly on the exact value of the mass associated to a 
given T e s and therefore these discrepancies are almost entirely a 
consequence of the different luminosities predicted by theoreti- 
cal model, for any particular T e g. 

To illustrate this point, we have computed synthetic photom- 
etry for the 1 Myr theoretical isochrones of DM98, and the 2 Myr 
isochrones of SDF00, which are adequate to represent the typical 
ONC population (see D10). The different ages of these models 
are likely due to different age scales, as discussed above. 

We also consid er the 2 Myr old DUSTY model 
dChabrier et alj|20"o"ol) of the Lyon group, well suited for cool 
atmospheres (T e ff <2900 K). In order to cover the spectral range 
of our selected sample of stars, we extend the DUSTY model 
to hotter temperatures using the BCAH98 model with mix- 
ing length parameter a=1.5, which, a posteriori, provides the 
best match between theoretical predictions and our empirical 
isochrone. 

We used the synthetic spectra of All ard et all (HoiO) and 
the empirical corrections to the atmospheric model derived in 
Sect. [4] In Fig. [13] we plot H, (J-H) and (H-K s ) as functions of 
r e ff for the various models. 

Figure [13] shows that there are significant differences in the 
H-T e g relation. The DUSTY curve provides the closest match to 
our empirical model down to spectral type M5; for later types, 
it drastically turns down such to overpredict logr e jf by ~0.4 
(corresponding to ~ 250 K). The same discrepancy is found by 
Dunu v et~a"D (1201 Ol) . who compared the DUSTY models with 
the NIR spectra of late-M dwarfs with high-quality dynamical 
masses. Based on their findings, they suggest that the atmo- 
spheric models are more likely than the evolutionary models to 
be the primary source of the discrepancy, since roughly the same 
r e ff offset is observed over a wide range of masses, ages, and ac- 
tivity levels but the same temperature range. 

Small differences are present between the DUSTY and 
SDF00 isochrones in the overlapping spectral range, with the 
noticeable exception that the latter shows a bump in luminos- 
ity for M5 stars. On the other hand, the DM98 model is slightly 
flatter: it overestimates the temperature of late K stars and un- 
derestimates that of mid-M types. Then, at very late spectral 
types, it steeply decreases reconciling with the BCHA98 model. 
This because the DM98 isochrones drastically turn down in lu- 
minosity for log7 ' e ff<34 5, as shown in Fig. [TJ As discussed by 
Iffillenbrand et al.1 (|2008), this drop is a common feature of the 
PMS evolutionary models, which tend to under-predict the in- 
trinsic luminosities of low-mass stars. 
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Figure [13] also shows that despite the discrepancies between 
the various H-T^ relations, the situation is much better in what 
concerns synthetic colors. In this case, model predictions are 
consistent within ~0.03 mag. In particular, the synthetic J-H col- 
ors of the DM98 and SDFOO models are redder by ~0.03 mag 
than our empirical isochrone. This is due to the fact that, for 
early to mid M types, the masses provided by these models 
are ~40% lower than the corresponding masses given by the 
DUSTY model, and this leads to the color shift discussed in 
Sect. 14.1.21 This confirms that the discrepancies are mainly due 
to different predictions of the stellar masses, which reflects in the 
predicted H magnitudes. 

Our empirical H, (J-H) and (H-Ks) are also shown in Fig. [13] 
as functions of T e g . The H-T^s curve provides higher fluxes for 
spectral types >M7 than the DUSTY and DM98 models. In gen- 
eral, no model provides H magnitudes nicely consistent with our 
empirical curve. On the other hand, our empirical colors are fully 
consistent, within the error bars, with the model predictions. 

6.2. Color-Magnitude Diagrams 

Figure [14] summarizes these relations in the CMDs, together 
with the density distribution of the extinction-corrected photom- 
etry of our sample of stars. The agreement between the various 
isochrones is generally poor, as the predicted magnitudes may 
differ by up to 1 mag. 

The synthetic colors corrected for the inaccuracies of the at- 
mospheric models (Sect. [4] solid lines in Fig. [T4l improve over 
the original theoretical predictions (dashed lines), tracing more 
closely the color distribution of the ONC sources in the CMDs. 
This correction, however, does not fix the systematic under- 
prediction of absolute magnitudes. Our empirical isochrone, de- 
rived from a subset of these data, also traces the color distri- 
butions of the sources in both CMDs and provides brighter H 
magnitudes for late type stars. 
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Fig. 15. Top panel - Mass-r e ff relation provided by the models 
of DM98, DUSTY and SDFOO. We also show our smooth rela- 
tion. Bottom p anel - Luminos i ty func tion derived from the mass 
distribution of Mu ench et al.l d2002l) . and assuming our semi- 
empirical mass-r e ff relation shown in the top panel. 



7. The accuracy of the photometrically-determined 
IMF of the ONC 

The discrepancies between the intrinsic magnitudes predicted by 
different evolutionary models affect the derivation of the main 
physical parameters of individual stars and clusters. In this sec- 
tion we analyze their impact on the reconstruction of the IMF 
of a young cluster like the ONC and the possible improvement 
attained by using our empirical magnitudes and colors. To this 
purpose, we run an artificial experiment assuming a mass distri- 
bution of the ONC and our empirical isochrone to produce "true" 
NIR magnitudes. Then, we recover the luminosity function (LF) 
of the cluster and the IMF by using the various model relations. 

Let us remind first that in order to analyze the IMF of a clus- 
ter like the ONC (for simplicity, we shall assume that the mass 
distribution of a young cluster still reflects the IMF), a com- 
mon approach is to collect NIR photometric data to derive, after 
some assumptions on the age, distance and reddening, the in- 
trinsic LF jL of the population in a generic photometric band m. 
The next step is to convert the LF of the sample into the under- 
lying mass distribution using an appropriate Mass-Luminosit 
Relation (MLR), m (M) (see e.g. iD'Antona & Mazzitellilll99 
lLada & La da 2003). The mass distribution -rn can then be de 
rived as 

dN dN dm 

— (AO = -r (m M)) • — M). (2) 
dM dm dM 

In most cases, the MLR is based on theoretical predictions. 



In order to simulate a statistically significant stellar sam- 
ple, we randomly generate 10 5 stars with mass distributed ac- 
cording to the mult i-part power-law mass function derived by 
iMuench et all d2002h for the ONC (Fig. [17}, and extending the 
trend at very low masses down to M ~ 0.02M Q (thus neglecting 
the secondary peak they report at sub-stellar masses). 

We then convert masses into effective temperatures, and then 
into photospheric / magnitudes. We focus on the /-magnitudes 
because this band is usually assumed to be the most suitable 
one (together with the I and H bands) in deriving the IMF of 
young clusters, being relatively unaffected by ac cretion, disk 
emission and other circums tellar activity (see e.g. iHillenbrandl 
l!997blLucas & Rochell2000l) . 

Our empirical isochrone performs the second step of the pro- 
cess, converting T e ff into magnitudes, but does not provide any 
information on the first passage, the mass vs T e g relation. To this 
purpose, we derive an average mass-r e ff relation by combining 
the relations of the different theoretical models (Fig. [15}. 

Combining the mass-r e ff relation with the empirical 
isochrone listed in Table [2] we obtain our semi-empirical MLR 
(Fig. [Toll needed in Eq.[2] We extend it to spectral types earlier 
than K6 using the 1 Myr old model of DM98, which smoothly 
connects to our isochrone. The J-LF of our artificial sample of 
stars is shown in Fig. [TBI 

The last passage of our exercise requires converting the LF 
in Fig.[l5]back into a mass function by means of Eq.[2l assum- 
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Fig. 13. H, (J- H) and (H-Ks) as fu nctions of T e ff (from top to bottom). At the top of each panel we report the spectral type-Jeff 
conversion of iLuhman e t al. (2003]). In each panel, dots represent the spectral types listed in Table [2] for which the correction is 
available (see Table[T|i. The thin gray solid lines represent the uncertainties on our derived empirical colors, as shown in Fig.fTUl 



ing different MLRfl As a sanity check, we start with our semi- 
empirical MLR (Fig. [15}, used to generate the LF. Figure [17] 
shows that, besides some numerical oscillation, we recover our 
assumed mass distribution. This proves the robustness of our al- 
gorithm. We then perform the same computation assuming the 
/-band MLRs provided by the models discussed in Sect. [6] We 
apply a kernel smoother to the MLRs in order to reduce the 
numerical fluctuations induced by the limited and uneven mass 
samplings of the model relations (Fig.fToTl, obtaining the results 
presented in Fig.fTTI 

Let's focus first on the stellar regime. Using the theoretical 
models, we generally recover the main features of the simulated 
IMF: it increases with decreasing masses for M > O.6M and 
then flattens and reaches the maximum frequency. The recovered 
peak location depends on the assumed model: the DM98 and 
DUSTY models return the same characteristic mass (M = 0. 12 Q ) 
as the simulated IMF, while the model of SDF00 returns a 
slightly larger mass (M ~ O.23M ) corresponding to the peak 
position. 
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2 We remind the reader that the luminosity spread discussed in 
Sect. [5] prevents any luminosity-mass conversion on a star-by-star ba- 
sis. Nonetheless, this statistical approach remains valid when using the 
most appropriate MLR on a large sample of stars, as the errors should 
average down. 



Fig. 16. Smoothed MLRs (top panel) and the corresponding 
derivatives (bottom panel) for the theoretical models discussed 
in the text. 



Following Mue nch et all (120021) . we fit the IMFs using the 
power-law 

dN 



dM, 
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Fig. 14. (HJ-H) and (H,H-Ks ) CMDs (top and b ottom panel respecti vely) of the ONC, compared to the 1 Myr ol d isochrone mode l 

of DM98 (cyan line), the 2 Myr old isochrone of lChabrier et all (120001) (green line) and the 2 Myr old isochrone of lSiess et all d2 000) 
green line). Synthetic photometry for these models has been computed using the grid of atmospheres provided by lAllard et al.l 
2010 ) (in dashes) and then corrected using the corrections in Table[T](solid lines), assuming Ry=3.1. Dots mark the position of the 

spectral types listed in Table [2] along the isochrones (from K6 to M8.5 with increasing magnitude, from K6 to M6 for the model of 

SDF00). The dotted lines show the iso-density contours for the analyzed stellar sample. 



12 



G. Scandariato et al.: Empirical Near Infrared colors for low-mass stars and brown dwarfs in the Orion Nebula Cluster 



Table 3. Fit coefficients F derived assuming different model 
MLR 





Simulated 


DM98 


DUSTY 


SDF00 


M > O.6M 


-0.78" 


-1.07 


-2.19 


-0.77 


0.12M G < M < 0.6M o 


-0.15" 


-0.26 


-1.11 


-0.25 



Notes. (a) From lMuench et al. 12002). 



comparison it is evident that at the low-mass end, the IMF 
derivations appear to diverge, and it is unclear how to make de- 
tailed comparisons between the different methods. In any case, 
since photometrically-derived IMFs primarily depend on the the- 
oretical MLR extracted from the PMS model, the differences be- 
tween them likely stem from the inaccuracies in the evolutionary 
models at subsolar masses, as we show in Fig [171 



0.5 




0.5 0.0 -0.5 -1.0 -1.5 

Log(M/M Q ) 

Fig. 17. Comparison of the IMFs derived using different MLRs. 
The solid black line is the mass distribution of our artificial sam- 
ple. The dashed vertical line marks the separation between stars 
and brown dwarfs (M - O.O75M ). 

with turn-off mass at M — Q.6M G . In Table[3]we list the derived 
T coefficients: we find that using the theoretical MLRs we over- 
estimate the simulated F coefficient in both mass ranges. In par- 
ticular, the DUSTY model shows the largest deviation from the 
expected results. This is due to the fact that while the derivatives 
of the MLR provided by the other models are roughly constant 
in the M < 0.12M o , the derivative corresponding to the DUSTY 
model increases with decreasing mass. 

Considering the substellar domain, we find that both the 
DUSTY and the DM98 models provide an over-population of 
brown dwarfs, both of them generating a feature in the IMF (as 
th e one seen in the ONC I MF provided bv lMuench et al.1 d2002l) 
or iLucas & Rochd (2000)) which could be mistaken as a sec- 
ondary peak in the substellar regime. 

Splitting the IMFs in two parts around the peak (M = 
0.1 2M ) and integrating in the corresponding mass ranges, we 
find that the fraction of very low-mass objects (M < 0.1 2M ) 
compared to the numberof stars with O.12M < M < 1M is 1.4 
and 0.7, assuming the DM98 and DUSTY models respectively. 
These fractions are inconsistent with the true population of stars, 
which provides a 0.4 fraction. This because the two evolutionary 
models provide luminosities which decrease too rapidly with the 
stellar mass (see Fig. [TJ and, by consequence, the derivatives 
of the MLRs in the substellar regime are higher than our semi- 
empirical expectation. 

For the sake of comparison with observatio nal results, we 
remark that our Fig. [17] is similar to Fig. 18 in iMuench et alJ 
(12001 . where the authors summarize a few derivations of the 
IMF for the central 5'surrounding the Trapezium stars. Such 
a small area is projected against the thickest part of the back- 
groun d Orion Molecular Cloud (Ay > 25, Scandariato et alJ 
201 1), thus galactic background contamination is negligible. As 
a consequence, differences between the observationally-derived 
IMFs are due to differences between techniques (photometry, 
spectroscopy. . . ) and assumed theoretical models. From this 



8. Conclusion 

In this paper we have accurately analyzed a sample of ~300 stars 
with measured temperatures, luminosities and photospheric NIR 
photometry as a benchmark for current atmospheric and evolu- 
tionary models for low-mass PMS stars and brown dwarfs. 

We have compared the extinction-corrected photometry to 
the expected photo metry provided by the template spectra of 
lAllard et~ai] y010), finding that major improvements have been 
done in the synthesis of theoretical spectra. Nonetheless, we ob- 
tain indications on the lack of opacity in the H-band, likely due to 
the improper treatment of either the water vapor absorption pro- 
file or the collision induced absorption from H2. We thus propose 
the set of empirical corrections listed in Table [TJ to be regarded 
as additive terms to the synthetic NIR colors. These corrections 
are weakly influenced by the log g assumed to derive synthetic 
colors. 

We also analyzed the same sample of stars in order to de- 
rive the average isochrone of the ONC, reported in Table [2] The 
analyzed sample of stars show a magnitude spread of the order 
of ~3.5 mag, consistent the the ~1.5 dex luminosity spread re- 
ported in previous study of the ONC. This spread does not allow 
us to constrain the magnitude scale of the ONC isochrone and, 
by consequence, does not allow any mass determination based 
solely on the magnitude of stars. On the other hand, the J-H and 
H-Ks colors weakly depend on the photospheric luminosity, and 
are well constrained by our statistical analysis. 

Comparing our empirical isochrone to current theoretical 
models, we find that there is generally good agreement with the 
2 Myr old model of DUSTY and, to some extent, the 2 Myr old 
model of SDF00, both in magnitude (down to M5 types) and 
color scale. On the other hand, we find indications that the PMS 
isochrones of DM98 provide r e fj-log L relations slightly flatter 
than our observational data. 

Finally, we investigated how the theoretical models affect the 
photometric derivation of the IMF in the NIR domain. We find 
that PMS evolutionary models generally underestimate the in- 
trinsic luminosity of VLMSs and BDs, and this may lead to ar- 
tificial structures in the low-mass tail of young clusters' mass 
distribution. 

The empirical NIR colors we have derived, in 
the 2MASS system, can be readily converted to 
other photometric systems using the prescriptions at 
http : //www . ipac . caltech. edu/2mass/releases/allsky/doc/s 
Since low-mass stars slowly evolve in the HR diagram at very 
early stages, our empirical colors can conveniently be used as 
the intrinsic colors of young (age < 5 Myr) stellar clusters to 
derive NIR excess and extinction of individual stars. 
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